Diffuse Galactic Gamma Rays from Shock- Accelerated Cosmic Rays 
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A shock-accelerated particle flux oc p~ s , where p is the particle momentum, follows from simple 
theoretical considerations of cosmic-ray acceleration at nonrelativistic shocks followed by rigidity- 
dependent escape into the Galactic halo. A flux of shock-accelerated cosmic-ray protons with s ~ 2.8 
provides an adequate fit to the Fermi-LAT 7-ray emission spectra of high-latitude and molecular 
cloud gas when uncertainties in nuclear production models are considered. A break in the spectrum 
of cosmic-ray protons claimed by Neronov, Semikoz, & Taylor (PRL, 108, 051105, 2012) when 
fitting the 7-ray spectra of high-latitude molecular clouds is a consequence of using a cosmic-ray 
proton flux described by a power law in kinetic energy. 



Introduction. — One hundred years after the discovery 
of cosmic rays by Victor Hess in 1912 the sources of 
the cosmic radiation are still not conclusively established. 
Theoretical arguments developed to explain extensive ob- 
servations of cosmic rays and Galactic radiations favor 
the hypothesis that cosmic-ray acceleration takes place 
at supernova remnant (SNR) shocks [2j. A crucial pre- 
diction that follows from this hypothesis is that cosmic- 
ray sources will glow in the light of 7 rays made by the 
decay of neutral pions created as secondaries in collisions 
between cosmic-ray protons and ions with ambient mat- 
ter. The p + p — > 7T° — > 27 spectrum formed by isotropic 
cosmic rays interacting with particles at rest is hard and 
symmetric about e 7 = = 67.5 MeV in a log-log 

representation of photon number spectrum vs. energy Q. 

The AGILE and Fermi 7-ray telescopes have recently 
provided preliminary evidence for a 7r° — ► 27 feature 
in the spectra of the W44, W51C, and IC 443 SNRs 
4J. Whether this solves the cosmic-ray origin problem 
depends on disentangling leptonic and hadronic emis- 
sion signatures, including multi-zone effects [5j, and find- 
ing out if the nonthermal particles found in middle-aged 
SNRs as inferred from their spectral maps have the prop- 
erties expected from the sources of the cosmic rays. 

Supporting evidence that cosmic rays are accelerated 
at shocks comes from the diffuse Galactic 7-ray emis- 
sion. High-latitude Galactic gas separate from 7-ray 
point sources, dust and molecular gas provides a "clean" 
target for 7-ray production from cosmic-ray interactions. 
Because the cosmic-ray electron and positron fluxes are 
> 20 times smaller than the cosmic-ray proton flux at 
GeV energies [f|, the 7r°-decay 7-ray flux strongly dom- 
inates the electron bremsstrahlung and Compton fluxes. 
The diffuse Galactic 7-ray spectrum can be deconvolved 
to give the cosmic-ray proton spectrum, given accurate 
nuclear 7-ray production physics. Data from the Fermi 
Large Area Telescope (LAT) offer an opportunity to de- 
rive the interstellar cosmic-ray spectra unaffected by So- 
lar modulation 0, H[ . 

This problem is revisited in order to address a recent 
claim [8[ that fits to the diffuse Galactic 7-ray emission 
of Gould belt clouds imply a break at T^br — GeV 
in the cosmic-ray proton spectrum that represents a new 



energy scale. Note that a break in the kinetic-energy 
representation of the proton spectrum has been reported 
before @ . In this Letter we show that when uncertainties 
in nuclear production models are taken into account, a 
power-law momentum spectrum favored by cosmic-ray 
acceleration theory provides an acceptable fit to Fermi- 
LAT 7-ray data of diffuse Galactic gas, and produces a 
break in a kinetic energy representation at a few GeV 
from elementary kinematics. Comparison of theoretical 
7-ray spectra with data supports a nonrelativistic shock 
origin of the cosmic radiation, consistent with the SNR 
hypothesis. 

Production spectrum of 7-rays from cosmic-ray 
collisions. — The 7-ray production spectrum divided by 
the hydrogen density uh is given, in units of (s-GeV) -1 , 
by 
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Here and below, j(T p ) is the cosmic-ray proton intensity 
in units of cosmic-ray protons (cm 2 -s-sr-GeV) _1 , and k, 
the nuclear enhancement factor, corrects for the compo- 
sition of nuclei heavier than hydrogen in the cosmic rays 
and target gas [l(|. The term dapp^^o ^, 2 -y{T p , e 7 )/de 7 is 
the differential cross section for the production of a pho- 
ton with energy e 7 by a proton with kinetic energy T p , 
in GeV, and momentum p in GeV/c. 

The much studied and favored model for cosmic-ray 
acceleration is the first-order Fermi mechanism, which 
was proposed in the late 1970s [ll|. Test particles gain 
energies by diffusing back and forth across a shock front 
while convecting downstream. The downstream steady- 
state distribution function is given by f(p) oc p~ 3r /( r_1 ), 
where p is the momentum and r is the compression ra- 
tio. Consequently, the particle momentum spectrum 
oc P 2 f{p) oc P~ A > where A = (2 + r)/(r — 1) is the 
well-known test-particle spectral index that approaches 2 
(equal energy per decade) in the limit of a strong nonrel- 
ativistic shock with r — ► 4 [12]. After injection into the 
interstellar medium with A = 2.1 - 2.2, characteristic 
of SNR shocks, cosmic-ray protons and ions are trans- 
ported from the galactic disk into the halo by rigidity- 
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FIG. 1: Production spectra of secondary 7-rays made in 
p-p collisions from the models of Dermer [15[ and Kamae et 
al. pjj, using the empirical demodulated cosmic-ray proton 
spectrum given by eq. ([3]). 



dependent escape [13|, which softens their spectrum by 
S ks 0.5 - 0.6 units, leaving a steady-state cosmic-ray 
spectrum dN/dp oc p~ s , where s = A + 5 = 2.7 - 2.8. 
Thus we consider a cosmic-ray flux 

j sh (T p ) ex P(T P )(^r) oc P{T P )\$L\ (^)«p(T,)-' . 

" "" '"' (2) 

Gamma-ray production from cosmic-ray /matter 
interactions. — Secondary nuclear production in proton- 
proton collisions is described by isobar formation at 
energies near threshold T p = 0.28 GeV and scaling 
models at high-energies T p 3> m p . Uncertainties remain 
in the production spectra at T p « few GeV, where most 
of the secondary 7-ray energy is made in cosmic-ray 
collisions. 

Fig. 1 compares two models for 7-ray production from 
p-p collisions using an empirical fit to the demodulated 
cosmic-ray flux meausred |14| in interstellar space given 
by 



jdem(T p ) = 2.2(T P + m p ) 



-2.75 



(3) 



[l5| . The model of Dermer [lj| describes low-energy res- 
onance production in terms of the A(1232) resonance 
through Stecker's model [3|, and high-energy production 
by the scaling model of Stephens and Badhwar 16], with 
the two regimes linearly connected between 3 and 7 GeV. 
The model of Kamae et al. [l?} is a parametric represen- 
tation of simulation programs, and includes contributions 
from the A(1232) isobar and V(1600) resonance cluster, 



> 

CD 

^ 10" ; 
>- 



10"' 



Shock-acceleration spectrum, s = 2.75 




• Abdo et al. (2009) 

Dermer (1986), k = 1.45 

Kamae et al. (2006), k = 1.45 

Dermer (1986), k = 1.84 

Kamae et al. (2006), k= 1.84 



0.1 



10 



e (GeV) 

Y 



FIG. 2: Fits to the Fermi-LAT spectrum of the differential 
7-ray emissivity of local neutral gas 0], employing a shock- 
acceleration spectrum, Eq. Q, with s = 2.75, for the cosmic- 
ray proton spectrum, and the models of Dermer [l5l | and Ka- 
mae et al. [ItJ for 7-ray production. The nuclear enhancement 
factor k takes the value of 1.45 and 1.84, as labeled. 



non-scaling effects, scaling violations, and diffractive pro- 
cesses. This model is represented by functional forms 
that are convenient for astrophysical calculations. 

Fig. [T] shows that the two models agree within 20% at 
e 7 > 100 MeV, but display larger differences at e 7 < 10 
MeV, where the 7-ray production is energetically insignif- 
icant. Both models are in general agreement in the high- 
energy asymptotic regime, with the Kamae et al. model 
giving fluxes larger by ps 13% due to the inclusion of 
diffractive processes. The spectral peak in a vF v rep- 
resentation occurs at 500 MeV - 1 GeV for this pro- 
ton spectrum. The most significant discrepancy in the 
two models is by ~ 30% at the pion production peak 
near 67.5 MeV 



18] 



Because of the different approaches 
of the models, and little improvement in nuclear data 
bases from laboratory studies between model develop- 
ment, this comparison indicates that our knowledge of 
the 7-ray production spectrum in p-p collisions is uncer- 
tain, at worst, by 30% near the pion-production peak, 
and is better than 15% at e 7 > 200 MeV. 

Fits to Fermi-LAT j-ray data — The Fermi-LAT data 
of the diffuse galactic 7 radiation are fit using a shock 
spectrum given by Eq. ©. Guided by the high-energy 
asymptote of Eq. ([3]) , we let 



jcn(T p ) = 2.2p~ 



(4) 



with s = 2.75. The data shown in Fig. [5] are from regions 
in the third quadrant, with Galactic longitude from 200° 
to 260°, and galactic latitudes |6| ranging from 22° to 
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FIG. 3: Fermi-LAT spectrum of the differential 7-ray 
emissivity of Gould belt clouds Q. It is fit with a shock- 
acceleration model for the cosmic-ray pro ton spectrum with 
s = 2.85, using the models of Dermer [15|] and Kamae et al. 
[I?! for 7-ray production. The models are normalized to the 
flux at 1 GeV. 



60°. There are no molecular clouds in these regions, the 
ionized hydrogen column density is small with respect to 
the column density of neutral hydrogen, and 7-ray emis- 
sion from point sources is removed. The linear increase 
of the emissivity as traced by 21 cm line observations 
allows residual galactic and extragalactic 7 radiation to 
be subtracted, leaving only the 7-ray emission resulting 
collisions of cosmic rays with neutral gas, allowing for an 
absolute normalization to be derived for the emissivity. 
As can be seen, the shock-acceleration spectrum gives an 
acceptable fit to the data, and restricts the value of k to 
be < 1.8. Given that there must be some residual cosmic- 
ray electron-bremsstrahlung and Compton radiations at 
these energies, the restriction on k could be larger. 

The spectrum in Fig.[3] from the analysis of Fermi- 
LAT data given in Ref. Q shows the averaged 7-ray flux 
from molecular clouds in the Gould belt. The derived 
flux, which extends to « 100 GeV, matches the spec- 
trum of diffuse Galactic gas emission used by the Fermi 
team in their analysis of the extragalactic diffuse 7-ray 
intensity [T^ | , supporting the assumption that the clouds 
are porous to cosmic rays. The cosmic-ray proton shock 
acceleration spectrum given by Eq. ((4]) with s = 2.85 is 
seen to give an adequate fit to the data in Fig. [3J given 
the nuclear physics uncertainties in 7-ray production. 

This cosmic-ray proton flux is compared in Fig. [4] with 
a shocked spectrum given by Eq. (QJ with s= 2.85, with 
two power-law kinetic energy spectra (one multiplied by 
/3), and the spectrum, eq. assumed to represent the 
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FIG. 4: Theoretical cosmic-ray proton fluxes compared with 
recent PAMELA, Fermi-LAT, ATIC-2, and HESS measure- 
ments of cosmic-ray p, He, C, Fe, and electron and positron 
fluxes [53]. The shock acceleration spectra are given by Eq. 
(Q with s = 2.75 and s — 2.85, a power-law kinetic energy 
flux j oc T~ 2 - 85 , and a flux j oc /3T~ 2 ' 65 are shown, in ad- 
dition to the demodulated cosmic-ray proton spectrum, eq. 
Q, and the best fit spectrum of Ref. [8J. Solar modulation 
accounts for the difference between the measured and theo- 
retical cosmic-ray proton fluxes at T < 10 GeV/nuc. Papers 
cited in [2(| report systematic errors for different experiments. 



demodulated local cosmic-ray proton spectrum. The best 
fit of the functional form used by Neronov et al. (2012) Q 
is also plotted, and is in accord with the shock spectrum, 
Eq. (|4|), when allowance is made for the large uncertain- 
ties in derived parameters. Indeed, the allowed spectrum 
would be further limited by the Fermi-LAT data of Fig. 
[2] unless k < 1.2. The reduced \ 2 to determine goodness 
of fit should take into account uncertainties in nuclear 
production, and firm conclusions about a break depend 
on improved cross sections. 

Summary. — Fig. 4 compares theoretical and empirical 
local interstellar spectrum of cosmic-ray protons with re- 
cent measurements of the fluxes of cosmic-ray protons, 
He, C, Fe and electrons. The small fluxes of the heavier 
cosmic rays and electrons imply that 7-ray emission from 
these channels make a minor, but non-negligible contri- 
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bution to the 7-ray flux. 

The favored model for the origin of cosmic rays is non- 
relativistic shock acceleration by SNRs in the Galaxy. 
The simplest possible shock-acceleration model with pro- 
ton flux oc p~ s gives an adequate fit to Fermi-LAT data 
of high Galactic latitude gas and clouds, and is in ac- 
cord with expectations of a SNR origin for the cosmic 
rays. The use of a shock-spectrum for the cosmic-ray 
flux allows a range of calculations to be made that de- 
pend on knowing the low-energy cosmic-ray spectrum 
[2l| . The combination of a power-law momentum injec- 
tion spectrum and the observationally similar local in- 
terstellar spectrum potentially puts constraints on prop- 
agation models that have a large effect on the primary 
spectra, such as those involving strong reacceleration |6j . 

If cosmic rays are indeed accelerated by SNR shocks, 
then the ir° — >• 2-f feature from shock-accelerated cos- 
mic rays should be observed in the spectra of individ- 



ual SNRs. Indications for such a feature are found in 
a few middle-aged SNRs but final confirmation of 
cosmic-ray sources will require detailed spectral calcula- 
tions involving both shock-accelerated protons and lep- 
tons [22j], including radiative losses and escape, and com- 
parison with improving Fermi-hAT data resulting from 
increasing exposure and development of better analysis 
tools. 



I thank A. Atoyan, J. D. Finke, J. Hewitt, R. J. Murphy, 
and S. Razzaque for discussions, and Professor T. Kamae 
for supplying his secondary nuclear production code. I 
would also like to thank A. Strong for many interesting 
and useful comments, and I. Moskalenko and T. Porter 
for pertinent remarks. This work is supported by the 
Office of Naval Research and NASA through the Fermi 
Guest Investigator program. 



[1] For a historical review of the discovery of cosmic rays, see 
P. Carlson and A. de Angelis, European Physical Journal 
H, 35, 309 (2011). 

[2] V. L. Ginzburg and S. I. Syrovatskii, The Origin of Cos- 
mic Rays (New York: Macmillan, 1964). S. Hayakawa, 
Cosmic Ray Physics (New York: Wiley-Interscience, 
1969) (New York: Wiley-Interscience, 1969). V. S. 
Berezinskii, S. V. Bulanov, V. A. Dogiel, and V. S. 
Ptuskin, Astrophysics of Cosmic Rays (Amsterdam: 
North-Holland, 1990), edited by V.L. Ginzburg. 

[3] F. W. Stecker, Cosmic Gamma Rays (Baltimore, MD: 
Mono Book Co., 1971). 

[4] AGILE results for W44 are reported by A. Giuliani et 
al, Astrophys. J. Lett. 742, L30 (2011), and Fermi LAT 
results for W44, W51C, and IC 443 are given in the pa- 
pers by A. A. Abdo, et al, Science, 327, 1103 (2010), 
A. A. Abdo et al, Astrophys. J. Lett. 706, LI (2009), 
and A. A. Abdo et al, Astrophys. J. 712, 459 (2010), 
respectively. 

[5] For the physical requirement of multi-zone models in 
SNR spectral fitting, and the resultant ambiguity to dis- 
criminate leptonic and hadronic models, see A. Atoyan 
and C. D. Dermer, Astrophys. J. Lett. 749, L26 (2012). 

[6] I. V. Moskalenko and A. W. Strong, Astrophys. J. 493, 
694 (1998); see also A. W. Strong, I. V. Moskalenko and 
O. Reimer, Astrophys. J. 613, 962 (2004). 

[7] A. A. Abdo et al, Astrophys. J. 703, 1249 (2009). 

[8] A. Neronov, D. V. Semikoz, and A. M. Taylor, Phys. Rev. 
Lett., 108, 051105 (2012). 

[9] A. W. Strong, I. V. Moskalenko, and O. Reimer, Astro- 
phys. J., 537, 763 (2000), and GALPROP models since 
then. The break in the kinetic energy spectrum was also 
found by T. Delahaye, A. Fiasson, M. Pohl, and P. Salati, 
Astron. Astrophys., 531, A37 (2011). 
[10] As discussed in [3, this factor ranges from 1.45 - 2.0, 
with a value of 1.84 favored in the work of M. Mori, 
Astrophys. J., 478, 225 (1997). 
[11] A. R. Bell, Monthly Not. Roy. Astron. Soc, 182, 147 
(1978); A. R. Bell, Monthly Not. Roy. Astron. Soc, 
182, 443 (1978); R. D. Blandford, and J. P. Ostriker, 



Astrophys. J. Lett. 221, L29 (1978); G. F. Krymskii, 
Akademiia Nauk SSSR Doklady, 234, 1306 (1977). 
[12] For reviews of nonrelativistic shock acceleration, see R. 
Blandford and D. Eichler, Phys. Rept., 154, 1 (1987), 
L. O. Drury, Rpts. Prog. Physics 46, 973 (1983), and 
J. G. Kirk, 1994, Saas-Fee Advanced Course: Plasma 
Astrophysics, 24, 225 (New York: Springer, 1994). Ap- 
plications to 7-ray studies of SNRs are reviewed by S. P. 
Reynolds, Ann. Rev. Astron. Astrophys. 46, 89 (2008) 
and D. Caprioli, J. Cosmology Astropar. Phys. 5, 26 
(2011). 

[13] Cosmic-ray convective and diffusive escape from the 
Galaxy in terms of rigidity (momentum per unit charge) 
to explain the ratio of isotopes in the cosmic-rays was 
considered, e.g., by F. C. Jones, Astrophys. J. 229, 747 
(1979). 

[14] J. A. Simpson, Ann. Rev. Nucl. Particle Sci. 33, 323 
(1983) 

[15] C. D. Dermer, Astron. Astrophys. 157, 223 (1986). This 
model is developed more fully in the papers by R. J. Mur- 
phy, C. D. Dermer, and R. Ramaty, Astrophys. J. Suppl., 
63, 721 (1987) and R. C. Berrington and C. D. Dermer, 
Astrophys. J. 594, 709 (2003). See also C. Pfrommer and 
T. A. Enfilin, Astron. Astrophys. 413, 17 (2004; erratum: 
C. Pfrommer and T. A. Enfilin, Astron. Astrophys. 426, 
777, 2004). 

[16] S. A. Stephens and G. D. Badhwar, Astrophys. Space 
Sci., 76, 213 (1981). 

[17] The parametrization of the model of Kamae et al. is given 
by T. Kamae, N. Karlsson, T. Mizuno, T. Abe, and T. 
Koi, Astrophys. J. 647, 692 (2006; erratum in T. Kamae 
et al., Astrophys. J. 662, 779, 2007), and is based on the 
treatment of diffractive effects and scaling violations in 
the paper by T. Kamae, T. Abe, and T. Koi, Astrophys. 
J. 620, 244 (2005). See also N. Karlsson and T. Kamae, 
Astrophys. J., 674, 278 (2008). 

[18] For other secondary production models, see I. Cholis, M. 
Tavakoli, C. Evoli, L. Maccione, and P. Ullio, J. Cos- 
mology Astropar. Phys., 5, 4 (2012) l|arXiv:1106.5073|l . 
using Pythia, C.-Y. Huang, S.-E. Park, M. Pohl, and 



5 



C. D. Daniels, Astropar. Physics, 27, 429 (2007), who 
use the Kamae et al. model at low energies, and H. Sato, 
T. Shibata, and R. Yamazaki, larXiv:1205.2i"45l (2012), 
who compare with LHC data. 

[19] A. A. Abdo et al, Phys. Rev. Lett. 104, 101101 (2010). 

[20] Ab09: A. A. Abdo, et al., Phys. Rev. Lett., 102, 181101 
(2009); AclO: M. Ackermann, et al., Phys. Rev. D, 82, 
092004 (2010); Adll: O. Adriani, et al., Phys. Rev. 
Lett., 106, 201101 (2011); Aril: O. Adriani, et al., 
Science, 332, 69 (2011 ); Pa06: A. D. Panov, et al., 
|arXiv:astro-ph/0612377l (2006); Ah08: F. Aharonian, et 
al., Phys. Rev. Lett., 101, 261104 (2008); Ah09: F. 
Aharonian, et al., Astron. Astophys., 508, 561 (2009). 
See http://www.mpe.mpg.de/~aws/propagate.html for 
cosmic-ray data and resources. 

[21] For positron production in SNRs, see Y. Ohira, K. Kohri, 
and N. Kawanaka, Monthly Not. Roy. Astron. Soc. 421, 
L102 (2012). For ionization of gas by low-energy cosmic 
rays, see B. B. Nath, N. Gupta, and P. L. Biermann, 



2012, arXiv: 1204.4239 , who also note that a momentum 
spectrum is theoretically preferred. See R. Ramaty, B. 
Kozlovsky, and R. E. Lingenfelter, Astrophys. J. Suppl. 
40, 487 (1979), for the Galactic 7-ray line spectrum in- 
duced by low-energy cosmic-ray interactions. 
[22] J. Fang and L. Zhang, Monthly Not. Roy. Astron. Soc. 
384, 1119 (2008); Y. Ohira, K. Murase, and R. Ya- 
mazaki, Monthly Not. Roy. Astron. Soc. 410, 1577 
(2011). 

Note added in proof: A recent paper by Blasi, Am- 
ato, and Serpico [Phys. Rev. Lett. 109, 061101 (2012)] 
explains a supposed low-energy break in a power-law in- 
jection momentum spectrum as a result of advective ef- 
fects on cosmic-ray propagation (compare Ref. [13]). M. 
Kachelriess and S. Ostapchenko (arXiv:1206.4705) claim 
to find deviations from a cosmic-ray power-law momen- 
tum spectrum at low energy, but their treatment is sus- 
ceptible to some of the same criticisms as presented here. 



